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Abstract 

The principal paradigm for gamma-ray bursts suggest that the prompt transient 
gamma-ray signal arises from multiple shocks internal to the relativistic expansion. 
This paper illustrates some properties of diffusive acceleration at relativistic shocks 
that pertain to GRB models, providing interpretation of the BATSE/EGRET data. 
Using a standard Monte Carlo simulation, computations of the spectral shape, and 
the range of spectral indices are presented, as functions of the shock speed, magnetic 
field obliquity and the type of particle scattering. It is clear that while quasi-parallel, 
relativistic shocks with fields approximately normal to the shock plane can efficiently 
accelerate particles, highly oblique and perpendicular ones cannot unless the particle 
diffusion is almost isotropic, i.e. extremely close to the Bohm limit. Accordingly, 
an array of distribution indices should be present in the burst population, as is 
exhibited by the EGRET data, and even mildly-relativistic internal shocks require 
strong field turbulence in order to model the > 100 MeV observations. In addition, 
recent spectral fitting to burst data in the BATSE band is discussed, providing 
probes of the efficiency of injection of non-thermal electrons. 
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1 Introduction 



Gamma-ray bursts are among the most interesting and exotic phenomena 
in astrophysics. In standard gamma-ray burst (GRB) models, the rapidly ex- 
panding fireball cools, converting the internal energy of the hot plasma into ki- 
netic energy of the beamed, relativistically moving ejecta and electron-positron 
pairs. At the point where the fireball becomes optically thin and the GRB we 
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see is emitted, the matter is too cool to emit gamma-rays unless some mecha- 
nism can efficiently re-convert the kinetic energy back into internal energy, i.e., 
unless some particle acceleration process takes place. Diffusive shock (Fermi) 
acceleration is widely believed to be this mechanism (e.g., Rees & Meszaros 
f992; Piran 1999). The large energy release in radiation from GRBs, coupled 
to the limits on available energy from likely sources such as supernovae or 
coalescing compact objects (see Piran 1999 and Meszaros 2002 for reviews), 
requires that this acceleration of particles at shocks be efficient. Moreover, the 
rapid variability of the prompt emission together with the impressive power 
of these sources strongly suggests that their environs are moving ultrarela- 
tivistically (e.g. Paczyhski 1986), an indication reinforced by inferences from 
77 transparency arguments (e.g. Baring & Harding 1997). Therefore, a de- 
tailed understanding of particle acceleration at relativistic shocks is clearly 
motivated for GRB studies. 

Acceleration of particles at non-relativistic shocks has been extensively in- 
vestigated in the contexts of supernova remnant and heliospheric shocks. Dif- 
fusive acceleration at relativistic shocks is less exhaustively studied than that 
for non-relativistic flows, yet it may occur in extreme objects such as pulsar 
winds, hot spots in radio galaxies, jets in active galactic nuclei and micro- 
quasars, and GRBs. Early work on relativistic shocks was mostly analytical in 
the test-particle approximation (e.g., Peacock 1981; Kirk & Schneider 1987a; 
Heavens & Drury 1988; Kirk & Webb 1988). Some analytic work (Schneider & 
Kirk 1987; Baring & Kirk 1991) has explored nonlinear, cosmic ray modified 
shocks. Complementary Monte Carlo techniques have been employed for rela- 
tivistic shocks by a number of authors, including test-particle analyses by Kirk 
& Schneider 1987b and Ellison, Jones & Reynolds (1990) for parallel, steady- 
state shocks, and extensions to include oblique magnetic fields by Ostrowski 
(1991), Ballard & Heavens (1992) and Bednarz & Ostrowski (1998). 

A key characteristic that distinguishes relativistic shocks from their non- 
relativistic counterparts is their inherent anisotropy due to rapid convection 
of particles through and away downstream of the shock, since particle speeds 
v are never much greater than the downstream flow speed u 2 ~ c/3 . Accord- 
ingly, the diffusion approximation, the starting point for virtually all analytic 
descriptions of acceleration at non-relativistic shocks, cannot be invoked since 
it requires nearly isotropic distribution functions. Hence analytic approaches 
prove more difficult for relativistic shocks, though advances in special cases 
such as the limit of extremely small angle scattering (pitch angle diffusion) are 
possible (Kirk k Schneider 1987a; Kirk, et al. 2000; Keshet & Waxman 2004). 
This paper explores some of the distinctive properties of particle acceleration 
at relativistic shocks that are germane to the gamma-ray burst paradigm, with 
a focus on spectral issues, specifically the slope of the power-law distribution 
and the efficiency of generation of this non-thermal component. 
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2 Acceleration at Relativistic Shocks in GRBs 

A most attractive feature of non-relativistic shock acceleration theory is 
that the distribution of accelerated particles is scale-independent, i.e. a power- 
law. This is a consequence of high energy particles (those with speeds ii>ui, 
with u\ being the upstream flow speed) attaining isotropy in all pertinent ref- 
erence frames. At such energies, the principal transport equation describing 
the acceleration process, the diffusion-convection equation, can be solved an- 
alytically for plane shocks (e.g., Blandford & Ostriker 1978; Jones & Ellison 
1991), yielding the well-known result for the momentum distribution 

-j- oc p with a = T . (1) 

dp r — 1 v ' 

where r = ui/u 2 is the shock (velocity) compression ratio, p is the momen- 
tum. Eq. (1) is a steady-state, test-particle result, and the index a depends 
only on the velocity compression ratio r , i.e. hydrodynamic quantities. This 
elegant result does not carry over to relativistic shocks because of their strong 
plasma anisotropy. As a consequence, while power-laws are in fact created, 
the index a becomes a function of the flow speed, the field obliquity, and the 
nature of the scattering, all of which intimately control the degree of particle 
anisotropy. 

In the specific case of parallel (i.e. those where the field is locally normal 
to the shock surface), ultrarelativistic shocks, the analytic work of Kirk et al. 

(2000) demonstrated that as r\ = 1/yl — (ui/c) 2 — > oo , the spectral index 
a asymptotically approached a constant, a — > 2.23 (see also Keshet & Wax- 
man 2004), a value realized once I\ <^ 10 . This captivating result has been 
confirmed by Monte Carlo simulations (Bednarz & Ostrowski 1998; Baring 
1999, Achterberg et al. 2001; Ellison & Double 2002), and corresponds to the 
special case of compression ratios of r = 3 and the particular assumption 
of small scattering (pitch angle diffusion), specifically for incremental changes 
# scatt in a particle's momentum with angle # scat t <C 1 /r\ . Here we explore 
departures from this particular case that are appropriate for burst studies. 

2. 1 Relativistic Shocks: an Array of Spectral Indices 

The spectral index of the power-law distribution is a declining function 
of the Lorentz factor r\ for a fixed compression ratio, a characteristic evi- 
dent in Kirk & Schneider (1987a), Ballard & Heavens (1991), and Kirk et al. 
(2000) for the case of pitch angle scattering, and a property that extends to 
large angle scattering (Ellison, Jones & Reynolds 1990; Baring, 1999; Ellison 
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& Double 2004). Faster shocks, if parallel, generate flatter distributions if r is 
held constant (e.g. see the parallel shock, Bn i — 0° case in Fig. 2), a conse- 
quence of the increased kinematic energization occurring at relativistic shocks. 
A tabulation of this property, namely a values for different r in the case of 
pitch angle diffusion, is given in Baring (2004), and analytic approximations 
are derived in Keshet & Waxman (2004). Note that imposing a specific equa- 
tion of state such as the Jiittner-Synge one renders r a function of Ti so 
that this monotonicity property can disappear, as evinced in Fig. 2 of Kirk 
et al. (2000; see also Keshet & Waxman 2004). Moreover, for significant field 
obliquities, the trend with shock speed is reversed, as is evident in Fig. 2. 

Astrophysical models usually invoke the canonical compression ratio r = 
3 , the well-known result for a relativistic, purely hydrodynamic shock possess- 
ing an ultrarelativistic equation of state. However, in cases where the magnetic 
field becomes dynamically important (e.g. the termination shock for the Crab 
pulsar wind: Kennel & Coroniti 1984), the strong fields can weaken magne- 
tohydro dynamic shocks considerably, just as they do in non-relativistic situa- 
tions. Moreover, in ultrarelativistic shocks, where pressure anisotropy can be 
significant, Double et al. (2004) observed that the shock could be strengthened 
or weakened, depending on the nature of the pressure anisotropy, which must 
be a significant function of the angle, Oeni , the upstream field makes to the 
shock normal. Hence, it is anticipated that a will be a function of Oeni • The 
influence of pressure anisotropy on the shock compression will be greatest in 
cases where the index is low, i.e. a ^ 2 , so that the accelerated particles can 
be influential in determining the dynamics of the (non-linear) shock. 

Also of interest is the fact that the slope of the nonthermal particle dis- 
tribution depends on the nature of the scattering, a feature evident in the 
works of Ellison, Jones & Reynolds (1990), Bednarz & Ostrowski (1998) and 
Baring (1999) . The asymptotic, ultrarelativistic index of 2.23 is realized only 
in the mathematical limit of pitch angle diffusion (PAD), where the particle 
momentum is stochastically deflected on arbitrarily small angular (and there- 
fore temporal) scales. In practice, PAD results when the scattering angle 6? scatt 
is smaller than the Lorentz cone angle 1 /Ti in the upstream region. In such 
cases, particles diffuse in the region upstream of the shock only until their 
angle to the shock normal exceeds around 1 /Ti . Then they are rapidly swept 
to the downstream side of the shock. 

Particle distributions obtained from the Monte Carlo simulation of acceler- 
ation at relativistic shocks developed by Ellison, Jones & Reynolds (1990) are 
exhibited in Fig. 1. They demonstrated that for large angle scattering (LAS, 
with # scat t ~ 7r ) the spectrum is highly structured and much flatter than E~ 2 . 
Such a case is exhibited in the Figure. The structure is kinematic in origin, 
where large angle deflections lead to distribution of fractional energy gains 
between unity and Y\ . Gains like this are kinematically analogous to the en- 
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Fig. 1. Particle distributions from a parallel ( 0Bni = 0) relativistic shock of 
r = 3 and Lorentz factor Ti = 5 , obtained from a Monte Carlo simulation (El- 
lison, Jones & Reynolds 1990; Baring, 1999). Scattering is modeled by randomly 
deflecting particle momenta by an angle # sca tt within a cone whose axis coincides 
with the momentum prior to scattering. Distributions are depicted for three cases, 
#scatt < 14° , corresponding to pitch angle diffusion (PAD), large angle scattering 
(LAS: # scat t < 7r » 1/Ti ), and an intermediate case (dotted histogram). The distri- 
butions are divided by the E~~ 2 - 23 ultra-relativistic power-law determined by Kirk 
et al. (2000). 



ergization of photons by relativistic electrons in inverse Compton scattering, 
and are much larger on average than those realized in PAD (see Gallant & 
Achterberg 1999; Baring 1999). Each structured spectral segment in Fig. 1 
corresponds to an increment in the number of shock crossings, successively 
from 1 — > 3 — > 5 — > 7 etc., as illustrated by Baring (1999), that eventu- 
ally smooth out to asympotically approach an index of a ~ 1.5 . Clearly, such 
highly-structured distributions have not been inferred from radiation emission 
in gamma-ray bursts or any other astrophysical objects. Note that the Fi — 5 
results depicted here are entirely representative of ultrarelativistic shocks. 
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An intermediate case is also depicted in Fig. 1, with # scat t ~ 4/r\ . The 
spectrum is smooth, like the PAD case, but the index is lower than 2.23. Mag- 
netic turbulence could easily be sufficient to effect scatterings on the order of 
these angles, a contention that becomes even more salient for ultrarelativis- 
tic shocks with Ti 3> 10 . Clearly a range of indices can be supported when 
#scatt is chosen to be of the order of 1/Ti , and the scattering corresponds to 
the transition between the PAD and LAS limits. Hence, it is expected that 
various astrophysical systems will encompass a range of scattering properties. 
Accordingly, the continuous and monotonically decreasing behavior of a with 
#scatt , as indicated in the exposition of Ellison & Double (2004), highlights the 
significant range of distribution indices in relativistic shocks. 

In relativistic shocks, the diversity of spectral indices is enhanced by trans- 
port properties orthogonal to the mean field direction. The diffusion of par- 
ticles across mean field lines, becomes a critical element in the discussion of 
oblique or perpendicular shocks. In non-relativistic shocks, when the upstream 
angle 0Bni of the field to the shock normal is significant, diffusion of parti- 
cles in the downstream region struggle to compete with convective losses, and 
transport back upstream of the shock layer becomes inefficient. This leads to 
the quenching of injection of thermal particles of speed v <: u x , which fail to 
return to the shock after one crossing to the downstream side. Accordingly, 
for @Bni ~ 30° , the Fermi acceleration process ceases when transport across 
field lines is suppressed (Baring, Ellison & Jones 1994), and for 0Bni < 30° , 
while power-law superthermal distributions are realized, their normalization 
is a strongly declining function of Bn i • 

In relativistic shocks, this phenomenon translates to a dramatic steepening 
of the distribution above thermal energies. When u\ ~ c , an oblique shock is 
inherently superluminal, so that convective losses are pervasive for all particle 
speeds, not just slightly suprathermal ones. Such losses must diminish the 
nonthermal population, and since the loss rate is purely a function of particle 
speed v (Peacock 1981; Jones & Ellison 1991), which is effectively pinned 
at v ~ c , and u<i , the overall effect is to increase the spectral index while 
retaining power-law character. Monte Carlo model indices are illustrated in 
Fig. 2, where the simulation output was acquired in the absence of cross field 
diffusion (i.e. for perpendicular spatial diffusion coefficient k± = \±v/3 = , 
corresponding to weak turbulence). Increasing ©Bni results in a rapid rise in 
a corresponding to a suppression of acceleration. Essentially, for ©Bni ^ 25° , 
acceleration is virtually non-existent for YiU\/c ^ 1 . Note that k\\ = A||i>/3 is 
the component of the spatial diffusion coefficient parallel to the field, so that 
defines the effective spatial scale of diffusion along the field lines. 

If the shocked plasma is strongly turbulent, the system will be driven 
towards the B ohm- diffusion limit, where diffusion coefficients are similar par- 
allel and perpendicular to the field, i.e. k± ~ k\\ , and transport is effectively 
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Fig. 2. Particle distribution indices a , for dn/dp oc , from oblique ( ©Bni > ) 
relativistic shocks of r = 3 as a function of T\(3i , obtained from a Monte Carlo 
simulation (Ellison, Jones & Reynolds 1990; Baring 1999) in the limit of pitch angle 
diffusion. Results are depicted for the case of zero diffusive transport perpendicular 
to the mean field, i.e., k± = for the component, perpendicular to B, of the spatial 
diffusion coefficient k = Xv/3 , where A is a particle's diffusive mean free path. The 
index is insensitive to G Bn i for non-relativistic shocks, but rapidly increases with 
obliquity for relativistic ones, underlining their inherent inefficiency. 



isotropic. Efficient transport across field lines returns the particles to the shock 
from the downstream region, effectively circumventing convective losses, and 
accordingly flattening the power-law distribution. In such cases, the strength 
of the turbulence renders quasi-perpendicular shocks much like parallel ones, 
and the distinction of field obliquity becomes less meaningful. It is anticipated 
that transport near the Bohm limit would be essential to generate a 3 , i.e. 
indices relevant for acceleration applications to burst models. This is borne 
out in Bednarz & Ostrowski (1998) and the recent work of Ellison & Double 
(2004), whose test-particle Monte Carlo simulation results of exhibited the 
expected monotonic decrease in a with an increase in k±/k\\ . 
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It must be remarked that the transport of particles in the shock environ- 
ment possesses a complex dependence on the nature of the dynamic fields. 
Niemiec & Ostrowski (2004) very recently explored the acceleration process 
at relativistic shocks by computing particle trajectories near a shock in the 
presence of injected, strong field turbulence. They found that (i) the spectral 
index was dependent on the strength of the turbulence, and more interestingly, 
(ii) that large perturbation amplitudes 5B/B <^ 1 could actually render a 
steepening of the distribution of accelerated particles for oblique, subluminal 
shocks. Niemiec & Ostrowski (2004) contended that this characteristic was 
largely a consequence of reduced reflection in the shock layer, and observed 
that the effect was generally reversed in oblique, superluminal shocks, where 
increasing SB/B flattened the distribution somewhat. Their investigations, 
and the Monte Carlo technique, of course do not explore quasi-coherent elec- 
trodynamic acceleration that results from charge separation potentials and 
currents in the shock layer; such plasma physics properties are best probed 
using particle-in-cell (PIC) simulations, which have been recently exploited by 
Silva et al. (2003) and Nishikawa et al. (2003) to study the development of the 
Wiebel instability and associated acceleration in relativistic plasma shocks. 
These injected turbulence and PIC simulations add to the evidence that the 
spectrum of particles accelerated at relativistic shocks is not universal. 



2.2 Modeling Observed GRB Prompt Emission 

The sensitivity of the power-law index to shock obliquity and the nature 
of turbulent transport immediately indicate that GRB spectra should pos- 
sess a diversity of indices, if a shock acceleration model is applicable to the 
burst paradigm. This is in fact manifested in the energy range well above the 
~ 1 MeV peak of emission in data taken from the EGRET experiment on 
the Compton Gamma- Ray Observatory (CGRO), where the half dozen or so 
bursts seen at high energies have a broad range of spectral indices (Dingus 
1995), namely a ~ 2 — 3.7 for dn/de^ oc e~ a . This result suffers from limited 
statistics due to (i) the nature of bursts, and (ii) to EGRET's field of view 
being more constrained than that for BATSE, the principal GRB experiment 
on CGRO. The GLAST mission will provide a more refined determination of 
the distribution of burst indices above 30 MeV after its launch in 2007. 

The EGRET data suggest that relatively flat indices are more common, 
though there is an obvious observational bias against distinguishing a ^ 6 
cases, since poor statistics at the uppermost energies will degrade index de- 
termination in such cases. If particle acceleration in bursts is indeed this effi- 
cient, then highly oblique shocks cannot be present, unless the turbulence is 
extremely strong; the results of Niemiec & Ostrowski (2004) support this con- 
tention. For prompt EGRET emission, the prevailing scenario is that shocks 
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internal to the GRB flow/blast wave are responsible for the dissipative con- 
version of bulk kinetic energy to observable radiation (e.g. see Piran 1999; 
Meszaros 2002). These shocks are necessarily mildly- relativistic, so a ~ 3 — 4 
values are reasonable for near the Bohm diffusion, even if the shocks are highly 
oblique. In contrast, optical afterglow models invoke the outer blast wave shock 
as the site for their energization, and in all probability, this ultrarelativistic 
shock is quasi-perpendicular. In such a shock environment, it is difficult to 
generate indices o ~ 3 even near the Bohm limit, so it is unclear how well 
Ti 3> 1 shocks will be able to model burst afterglow spectra. 

The power-law index is not the only acceleration characteristic germane 
to the GRB problem: the shapes of the particle distributions at thermal and 
slightly suprathermal energies are also pertinent. This energy domain samples 
particle injection or dissipational heating in the shock layer, and is readily 
probed for electrons by the spectrum of prompt GRB emission in the BATSE 
band. Tavani (1996) obtained impressive spectral fits to several bright BATSE 
bursts using a phenomenological electron distribution and the synchrotron 
emission mechanism. This has been a driver for the interpretation of burst 
spectra. There are issues with fitting low energy (i.e. ^ 100 keV) spectra in 
about 1/3 of bursts (e.g. Preece et al. 1998) in the synchrotron model, yet this 
emission mechanism still remains the most popular candidate today. Tavani's 
work did not directly address theoretical characteristics of distributions of 
accelerated particles at slightly suprathermal energies. 

Perspectives based on acceleration theory underpinned the recent analy- 
sis of Baring & Braby (2004), who pursued a program of spectral fitting of 
GRB emission using a linear combination of thermal and non-thermal electron 
populations. These fits demanded that the preponderance of electrons that are 
responsible for the prompt emission reside in an intrinsically non-thermal pop- 
ulation. This requirement strongly contrasts particle distributions obtained 
from acceleration simulations, exemplified by those depicted in Fig. 1. The 
consequence is obviously a potential conflict for acceleration scenarios where 
the non-thermal electrons are drawn directly from a thermal gas (the virtu- 
ally ubiquitous case), unless radiative efficiencies only become significant at 
highly superthermal energies. This does not necessarily mean that acceleration 
at relativistic shocks is not operating in bursts. It is possible that somehow, 
relativistic shocks can suppress thermalization of electrons, though such a con- 
jecture has no definitive simulational evidence to support it at present. A po- 
tential resolution to this dilemma is that strong radiative self- absorption could 
be acting, in which case the BATSE spectral probe is not actually sampling the 
thermal electrons. It is also possible that other radiation mechanisms such as 
pitch-angle synchrotron, or jitter radiation may prove more desirable. A goal 
of future work will be to ascertain whether a shock acceleration paradigm can 
be truly consistent with the GRB emission that is observed. 



9 



3 Conclusions 



This paper has discussed two key acceleration issues for models of gamma- 
ray burst sources, namely expectations for the power-law index a of the non- 
thermal population, and the efficiency for which these particles are injected 
into the acceleration process. The recent analysis of Kirk et al. (2000), ap- 
plicable for the specific case of pitch angle diffusion at plane parallel shocks, 
has demonstrated that the index of this power-law asymptotically approaches 
a = 2.23 as the Lorentz factor T\ of the upstream flow in the shock rest 
frame tends to infinity. Here it is illustrated how this widely-quoted result is 
not universal, and that an array of indices are possible in general. Results from 
a kinetic Monte Carlo simulation of diffusive shock acceleration are presented. 
These indicate that the value of a is sensitive to the nature of the scattering, 
becoming much smaller than 2.23 when the particle transport experiences 
angular deflections larger than around . The simulated distributions in 
I\ ^> 1 cases exhibit a rapid steepening (i.e., higher a) as the magnetic field 
obliquity Bn i of the shock rises above zero, in the particular case where 
diffusion perpendicular to the mean field is absent: efficient acceleration is 
effectively quenched in moderately oblique and quasi-perpendicular truly rel- 
ativists shocks unless there is strong cross field diffusion. The array of possible 
values of a is commensurate with EGRET data on a handful of bursts. This 
characteristic favors an internal shock model for the prompt emission, where 
the shocks are only mildly-relativistic in the comoving expansion frame, and 
the spread of a values is narrower. It is also noted that the acceleration 
simulations only inject a minority of thermal particles into the diffusive accel- 
eration process, contrasting electron distributions inferred in spectral fits for 
bright CGRO bursts. This poses a potential conflict for acceleration models 
that motivates further exploration of the energization of electrons via detailed 
hydrogenic and pair plasma simulations of shock acceleration. 
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